Finite-temperature chiral transitions in QCD with the Wilson quark 
action * 

S. Aokif Y. Iwasakif< b K. Kanayaf< b S. Kaya'r A. Ukawaf and T. Yoshie a < b 

institute of Physics, University of Tsukuba, Ibaraki 305, Japan 

b Center for Computational Physics, University of Tsukuba, Ibaraki 305, Japan 

institute of Particle and Nuclear Studies, High Energy Accelerator Research Organization(KEK), 
Tsukuba, Ibaraki 305, Japan 

We investigate the finite-temperature phase structure and the scaling of the chiral condensate in lattice QCD 
with two degenerate light quarks, using a renormalization group improved gauge action and the Wilson quark 
action. We obtain a phase diagram which is consistent with that from the parity-flavor breaking scenario. The 
scaling relation for the chiral condensate assuming the critical exponents and the scaling function of the three 
dimensional 0(4) model is remarkably satisfied for a wide range of parameters. This indicates that the chiral 
transition in two flavor QCD is of second order in the continuum limit. 



1. Introduction 

As a step toward the clarification of the finite- 
temperature QCD transition, it is important to 
investigate the nature of the transition on the lat- 
tice with two degenerate light quarks. In a previ- 
ous study using the Wilson quark action and a 
renormalization group (RG) improved gauge ac- 
tion 
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with Co = 1 — 8ci and c\ = —0.331 we studied 
the nature of the phase transition and the scaling 
behavior of the chiral condensate at (3 > (3 c t (~ 
1.35) on an N t = 4 lattice, where /3 c t is the value 
of P at the chiral transition point. In the present 
work, we extend the study to [3 < (3 c t performing 
simulations at j3 = 1.1, 1.2, and 1.35 on an 8 3 x 4 
lattice. We also determine the phase structure 
of the chiral limit, in particular, for [3 < [3 c t at 
N t — 4 on a (/3, K) plane. 
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Figure 1. Phase diagram for Np = 2 with the 
RG improved gauge action and the Wilson quark 
action. 

2. Phase diagram 

Fig.[l] shows our result for the phase diagram on 
a (/3,K) plane. In the previous workQ, the line 
of the zero-temperature chiral limit K C (T = 0) 
defined by the vanishing point of the pion mass 
was determined on an 8 4 lattice, and the finite- 
temperature transition/crossover line K t for Nt — 
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Figure 2. m£ and 2m g vs. 1/K for Nf = 2 ob- 
tained at (3 = 1.2 and 1.35 on an 8 3 x 4 lattice. 

4 was determined on an 8 3 x 4 lattice. The K t 
line crosses the K C (T = 0) line at f3 ct « 1.35. 

In this work, we estimate the location of the 
finite-temperature chiral limit K c (N t = 4) for 
N t = 4, defined by the vanishing point of the 
pion screening mass m„. In Fig.^, we plot rn^ 
and the quark mass m q for N t = 4 at fixed (3 
as a function of 1/K, where m q is defined by 
an axial Ward identity HQ. In the upper figure 
for (3 = 1.2, there are two chiral limits. On the 
other hand, at /3 = 1.35 ~ /3 c t in the lower fig- 
ure, we cannot find a clear gap between the two 
chiral limits. These results imply that the line 
K c (N t — 4) for the chiral limit turns back to- 
wards strong coupling around f3 c t, forming a cusp, 
as shown in Fig.|l|. This structure is consistent 
with that expected from the parity-flavor break- 
ing scenario^], which was confirmed for the stan- 
dard one plaquette gauge action and the Wilson 
quark action^]. 

3. Scaling analysis 

The magnetization M near the second order 
transition is expected to be described by a single 
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Figure 3. Subtracted chiral condensate (^ r ^ r ) su b 
as a function of 2m q a. 

scaling function; 

M/h 1 ' 5 = (2) 

with h the external magnetic field and t = [T — 
T c (h = 0)]/T c (h = 0) the reduced temperature. 

According to the universality argument, when 
the chiral transition is of second order, two fla- 
vor QCD belongs to the same universality class 
as the three dimensional 0(4) spin model0, and 
the chiral condensate should satisfy the scaling 
relation @ with the 0(4) exponents and the 
0(4) scaling function. We identify h = 2m q a, 
t = P-P c t, and M = (**) su b, where M is a 
subtracted chiral condensate defined through an 
axial Ward identity ||, 

(**) sub = 2m q a(2K f ^(tt^MO)). (3) 

X 

Our results for (^"I^sub are shown in FigJ^. 

We make a fit to the scaling function of the 
0(4) model by adjusting (3 ct and the scales for 
t and h, with the exponents fixed to the 0(4) 
values ||, including all the data in the range < 
2m q a < 0.8 and (3 < 2.0 shown in Fig.[| We 
note that the data for the lightest quark mass 
{m q = 0.06-0.07) at (3 = 1.1, 1.2, and 1.35 are 
slightly off the scaling curve. Therefore, we make 
a fit excluding these three points. The result of 
the fit with x 2 /df = 0.72 is shown in Figga). 
The result shows that scaling works well in the 
t < region as well as the t > region, except the 
three lightest quark mass data at t < 0, indicated 
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Figure 4. Fits to the scaling function with 
(a) 0(4) and (b) MF exponents. Solid curves 
are scaling functions obtained in an 0(4) spin 
modelQ and by a MF calculation, respectively. 

by the filled symbols. The adjusted (3 ct = 1.36(1) 
is consistent with f3 ct — 1.35(1) obtained by the 
previous fit using only the t > data|IJ. 

One possible origin of the deviation from the 
scaling curve of the lightest quark mass data at 
t < is a finite size effect, because, in the con- 
fining phase, we expect that finite size effects are 
severe at small quark mass. To explore this pos- 
sibility, we are performing simulations on L 3 x 4 
(L = 12 and 16) lattices for the lightest quark 
mass data at /3 = 1.1, 1.2 and 1.35. The statistics 
we have accumulated so far (about 50 configu- 
rations each) is not yet high enough to obtain a 
definite conclusion about the systematic size de- 
pendence of the deviation from the scaling curve. 

Another possible origin for the deviation is the 
explicit chiral breaking due to the Wilson term, 
which is expected to be large at small (3. It is 
plausible that this effect becomes visible as devi- 
ation from scaling, when the explicit chiral viola- 
tion due to the Wilson term becomes larger than 



that from the quark mass. The fact that we ob- 
serve the deviation from the scaling at the lightest 
quarks mass for small (3 region is consistent with 
this interpretation. To confirm this interpreta- 
tion it is necessary to make a simulation at larger 
Nt to reach a larger (3 region, or use an improved 
quark action. 

We also study the possibility of the MF scal- 
ing [^0). We perform the scaling test using MF 
exponents and the MF scaling function. Because 
m T at (9=1.5 does not vanish in the chiral limit 
(cf. Fig.3 in Ref.g), we restrict f3 ct < 1.5. We 
then obtain the best fit shown in Fig.^(b) with 
X 2 /df = 6.2, to be compared with x 2 /df = 0.72 
obtained for the 0(4) case [Fig.|(a)]. We exclude 
the MF scaling, as the data is much more scat- 
tered than in the 0(4) case. 

The success of the scaling with the 0(4) expo- 
nents, albeit with the three lightest quark data at 
t < excluded, indicates that the chiral transi- 
tion for two flavor QCD is of second order in the 
continuum limit. 

Numerical simulations are performed with Fu- 
jitsu VPP500/30 at the University of Tsukuba. 
This work is in part supported by the Grants- 
in-Aid of Ministry of Education (Nos. 09246206, 
09304029, 08640349, 08640350). 



REFERENCES 

1. Y. Iwasaki, K. Kanaya, S. Kaya, and T. 
Yoshie, Phys. Rev. Lett. 78 (1997) 179. 

2. Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
preprint of Tsukuba, UTHEP-118 (1983), un- 
published. 

3. M. Bochicchio et al, Nucl. Phys. B262 (1985) 
331. 

4. S. Itoh et al, Nucl. Phys. B274 (1986) 33. 

5. S. Aoki, Phys. Rev. D30 (1984) 2653; Phys. 
Rev. Lett. 57 (1986) 3136; Nucl. Phys. B314 
(1989) 79. 

6. S. Aoki, A. Ukawa and T. Umemura, Phys. 
Rev. Lett. 76 (1996) 873. 

7. R. Pisarski and F. Wilczek, Phys. Rev. 
D29(1984) 338; F. Wilczek, Int. J. Mod. Phys. 
A7(1992) 3911; K. Rajagopal and F. Wilczek, 
Nucl. Phys. B399(1993) 395. 



4 



8. K. Kanaya and S. Kaya, Phys. Rev. D51 
(1995) 2404. 

9. D. Toussaint, Phys. Rev. D55 (1997) 362. 

10. A. Kocic and J. Kogut, Nucl. Phys. B455 
(1995) 229. 



